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Abstract 

We give two "complementary" descriptions of the curve T parametrizing doublc- 
Bloch solutions to the difference analogue of the Lame equation. The curve depends 
on a positive integer number I and two continuous parameters: the "lattice spacing" 
rj and the modular parameter r. Apart from being a covering of the elliptic curve with 
the modular parameter r, F is a hyperelliptic curve of genus 11. We also point out 
connections between the spectral curve and representations of the Sklyanin algebra. 



1 Introduction 

Schrodinger operators with a periodic potential usually have infinitely many gaps in the 
spectrum. Exceptional cases, when there are only a finite number of gaps, are of particular 
interest for the theory of ordinary differential equations as well as for applications. Their 
study goes back to classical works of the last century. The renewed interest to the theory 
of finite-gap operators is due to their role in constructing quasi-periodic exact solutions to 
non-linear integrable equations. 

Among known examples of the finite-gap operators, the most familiar one is the classical 
Lame operator 

C = ~+£(£ + l)p{x), (1.1) 

where p(x) is the Weierstrass jp-function and I is a parameter. The finite gap property of 
higher Lame operators for integer values of I was established in jlj]. If i is a positive integer, 
then there exists a differential operator of order 2£ + 1 that commutes with C, so the Lame 
operator has exactly £ gaps in the spectrum. Such a remarkable spectral property is a 
signification of a hidden algebraic symmetry underlying the spectral problem. 

In [Q], a connection between the finite-gap integration theory of soliton equations and 
the representation theory of Sklyanin algebra was found and the following difference 
analogue of the Lame operator was proposed: 

_ 9 1 (2x-2£r ] ) ft (2s + 26?) dx 

e 1 (2x) + e 1 (2x) ' 1 j 
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Here 6i(x) = 9\{x\t) is the odd Jacobi ^-function, £ is a non- negative integer and rj G C is 
a parameter which is assumed to belong to the fundamental parallelogram with vertices 0, 
1, r, 1 + t. The origin of the operator ( |1.2| ) is traced back to Sklyanin's paper M of 1983, 
where a functional realization of the Sklyanin algebra was found. Namely, L coincides with 
one of the four generators of the Sklyanin algebra in the functional realization. Therefore, 
the Sklyanin algebra provides a natural algebraic framework for analyzing the spectral 
properties of the operator L. (A different algebraic approach to the difference analogues 
of the Lame operators was proposed in ||.) Nowdays, another face of this operator is 
probably more familiar: it is the hamiltonian of the elliptic two-body Ruijsenaars model 

!• 

As is already expected from the relation to the Sklyanin algebra, the spectral problem 
m> = E^fr is closely connected with the simplest one-site XYZ spin chain of spin £ at the 
site. Indeed, the operator L is proportional to the trace of the quantum L-operator of this 
model. Integrable spin chains of XYZ-type can be solved by the generalized algebraic 
Bethe ansatz , , || . In our case the Bethe ansatz approach amounts to looking for the 
eigenfunctions of the form 

£ 

=K x/r '\{e l (2x-2x j ), 

3=1 

where K and Xj are parameters. If they are constarined by the system of Bethe equations 

2 01(2^-217?) " giggj- 2 Xj -2r]) 
6 1 (2x i + 2£ V ) -2x^+2^' % ' 

then ^ is an eigenfunction of L. The energy E obtained from the eigenvalue equation at 
a particular value of x (say, x = £rj), 



e 1 (2£ V )f = \ 6 1 (2£r ] -2x J 



is then a multivalued function of K. It becomes single- valued on the spectral curve Y of the 
operator L, which, therefore, carries all the information about its spectral properties. The 
points P of the curve are solutions P = {K, xi, . . . ,xe} to the Bethe system. However, 
the description of the curve provided by the Bethe equations is neither the most economic 
nor very informative one (at least for small values of £). 

Let ^f(x) be a function on which the operator ( p..2[ ) acts. Putting \l/ n = ^>(nrj + Xq), 
we assign to (|1.2|) the difference Schrodinger operator L^/ n = A n ^/ n+ i + _B n \I/„_ 1 with 
quasiperiodic coefficients. The spectrum of a generic operator of this form has a structure 
of Cantor set type. If r] is a rational number, t] = P/(2Q), this operator has Q-periodic 
coefficients. In general, Q-periodic difference Schrodinger operators have Q stable bands 
in the spectrum. 

It was proved that for integer £ the operator L is algebraically integrable and, 
therefore, is a difference analogue of the classical Lame operator which can be obtained 
from L in the limit rj — > 0. Algebraic integrability of L implies, in particular, some 
extremely unusual spectral properties of this operator. Namely, the operator L given by 
eq. (|1.2|) for positive integer values of £ and arbitrary generic 77 has 2£ + 1 stable bands (and 
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2£ gaps) in the spectrum. Its Bloch functions are parametrized by points P = (w, E) of a 
hyperelliptic curve of genus 2£ defined by the equation 

w 2 = c J] (E 2 - El) , (1.3) 

i=i 

where c is a constant (introduced here for consistency with the definition of w given 
below). Moreover, eigenf unctions of the operator L at the edges of bands E = ±Ei span 
an invariant functional subspace for all generators of the Sklyanin algebra. 

The finite-gap property of the operator L means that there exists a difference operator 
W of finite order such that it can not be rerpesented as a polynomial function of L, 
and that commutes with L: [L, W] = 0. (This is the difference version of the Novikov 
equation for coefficients of the operators L and W.) It is well-known from the early days 
of finite-gap theory that the ring of operators commuting with the finite-gap operator is 
isomorphic to a ring of meromorphic functions on the corresponding spectral curve with 
poles at "infinite points". For difference operators this was proved in ||10|| , |TT|| . Therefore, 
the ring of operators commuting with L is generated by L and an operator W such that 

2£+l 

W 2 = c]l(L*- El) . (1.4) 

i=l 

The variable w in (|1.3|) is eigenvalue of the operator W, i.e., W^f = w^f, where \& is a 
common eigenfunction of L and W. The hyperelliptic curve ( |1.3| ) has two "infinite points" 
oo±, where the function E has first order poles. In for any algebraic curve with two 
punctures, a special basis in the ring of meromorphic functions with poles at the punctures 
was introduced. Due to the isomorphism between the ring of meromorphic functions with 
poles at oo ± and the ring of commuting operators, there exist operators commuting with 
L such that their eigenvalues on the common (Baker-Akhiezer) eigenfunction coincide 
with the basis functions. The explicit form of these operators and the operator W in 
the case when L is the difference Lame operator ( |1.2|) was found in [pl|. Remarkably, 



this commuting family (parametrized by a complex parameter) coincides with the Baxter 
Q-operator for the one-site XY Z -model with spin £. 

This paper is devoted to a more detailed analysis of the spectral curves of the difeffer- 
ence Lame operators for arbitrary positive integer values of £. Let us mention that spectral 
curves of the classical Lame operator (|1 . 1| ) and its Treibich-Verdier generalizations |L4] for 
small values of £ were studied in ||15||. Section 2 contains some algebraic preliminaries on 



quantum transfer matrices and the Q-operator for the simplest one-site XYZ-model with 
spin £ e Z + . This algebraic framework is very helpful since it allows one to represent the 
equation of the spectral curve in the most compact explicit form. The exposition in the 
first part of Section 3 follows that of the paper ||. Namely, we construct double-Bloch 
solutions to the difference Lame equation and obtain the spectral curve defined by two 
equations for three variables. One of these variables is the eigenvalue E, the other two 
parametrize the Bloch multipliers of the solution. Next, we show that E can be eliminated 
leaving us with one equation for two variables. In this realization, the fact that the spectral 
curve is a covering of the elliptic curve is transparent. However, the hyperelliptic property 
of the curve is implicit. In Section 4 we present the curve in the explicit hyperelliptic 
form. Finally, Section 5 contains a few results on the connection with representations of 
the Sklyanin algebra. 
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2 Quantum transfer matrices and the Q-operator for 
the one-site XYZ-model with spin £ 



This section contains selected ingredients of the quantum inverse scattering approach to 
XYZ spin chains Q specified to the case of the one-site "chain" with spin I. These 
constructions turn out to be particularly useful for analyzing the spectral curve of the 
difference Lame operator. 

We begin with a few formulas related to the Sklyanin algebra and its representations. 
Definitions and transformation properties of the Jacobi ^-functions 8 a (x\r) are listed in 
Appendix A. For brevity, we write 6 a (x\r) = 6 a (x). 

The elliptic quantum L-operator is the matrix 

/ 1 («)<s o + e 2 (u)s 1 + e 3 (u)s 2 \ 

L(«) = J (2.1) 

V e 2 (u)s l - e 3 (u)s 2 e 1 (u)s - e,(u)s 3 J 

with non-commutative matrix elements. Specifically, S a are difference operators in a com- 
plex variable x: 

_ e a+1 (2x-2£ V ) Qx _ 6 a+1 (-2x - 2£ V ) dx 
* a ~ 9 l {2x) 6 9i(2x) ' 1 j 

introduced by Sklyanin J| in 1983. Comparing with (|1.2|) , we identify L = So- The four 
operators S a obey the commutation relations of the Sklyanin algebra^ 

(— l) a+1 I a0 S a So = Ip^SpS-f — LfpS-fSp , 

(2-3) 

(— l) a+1 I a0 SoS a = IypSpS-f — Ip 7 S 7 Sp 

with the structure constants I a b = 9 a+ i(0)9b + i(2r]). Here a, b = 0, . . . , 3 and {a,/?, 7} 
stands for any cyclic permutation of {1,2,3}. The relations of the Sklyanin algebra are 
equivalent to the condition that the L-operator satisfies the "RLL = LLR" relation with 
the standard elliptic i?-matrix. 

The parameter i in ( p.2| ) is called spin of the representation. If necessary, we write 
S a = Sjp to indicate the dependence on i. When i £ |Z + , these operators have a finite- 
dimensional invariant subspace, namely, the space of even ^-functions of order 4£ 
(see Appendix A). This is the representation space of the (2£ + l)-dimensional irreducible 
representation (of series a)) of the Sklyanin algebra. 

Trace of L(w) (in the two-dimensional auxiliary space), that is the simplest quantum 
transfer matrix Ti(u), is proportional to L = Sq: 

T x {u) =trl_(¥) = 0i(u)<S o . (2.4) 



The whole family of commuting transfer matrices T s (u), s G Z + , is obtained from (|2.1| ) via 
the fusion procedure [0. We denote them by T s (u), s £ Z + . They commute for all values 

1 The standard generators of the Sklyanin algebra || are related to ours as follows: S a = 
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of s and u: [T s (u), T s /(u')] = 0. Here we do not need to recall the fusion procedure itself 
and refer the reader to [|17]] and [[0|], where integrable magnets of higher spin in the XX Z 
and XYZ case, respectively, were constructed by means of the fusion procedure. For our 
purpose it is enough to define the operators T s {u) by the recurrence relations (known as 
the fusion relations): 

s < 2£ : r 1 (w-sry)r s (M+ry)=T afl (ii)+0 1 (M-s7 ? -2^)0i(M-s77+2(£+l)77)T^ 1 (w+2ry), 
s = 2£ : Ti(u-2^)T m (u+7 ? ) = 0i(m)T 2W (u)+0 1 (u+2»7)0i(u-^)T 2< _i(w+2»7) , 
s > 11 : T 1 {u-srj)T s {u+r])=e 1 {u~sr]+2£ri)T^ 1 (u)+e 1 (u-sr]~2£ri)T s ^ 1 (u + 2r]) 

(2.5) 

with the "initial condition" Tq(u) = id, T\{u) = 9i(u)So. There is a useful determinant 
formula which solves the recurrence relations and represents T s (u) through Ti(u) [ |T8|j : 

T s {u) = det(T ii (s,u))i< i j< s , 0<s<2£, 

(s-U ^ 

T s (u) = \^j[9 1 (u + (2£+2i-s-l)r l )j det(T y (s, u))i<ij<. , s>2£, 
where 

T^(s,u) = 5 t ^ 1 9 1 (u + (s-2£-l-2i)ri) + 5 lJ+1 9 1 (u + (s + 2£+3-2i)ri) 

(2.7) 

+ (J ij -r 1 («+(s+l-2i)»7). 

Let \& be a common eigenfunction of L = iSo and T s (u), and let i£ be the eigenvalue of 
L: m> = Ety. Then the eigenvalue of T s (u) is a polynomial in i£ of degree s, which we 
denote by T s (u, E): 

\ T s (u)V = T a (u,E)V. 

The eigenvalues T s (u,E) are determined by eqs. fl2.6|), (|2.7D, where Ti(it+(s+l — 22)77) i s 
replaced by E9 1 (u+(s+l—2i)r]). As a function of it, T s (w, E) for 1 < s < 2£ is easily seen 
to belong to the space Q s of ^-functions of order s (for the precise definition see Appendix 
A). An important fact, not obvious from the definition ( |2.5|) , is that T s (u, E) for all s > 2£ 
belong to the space Q 2 e, i-e., the denominator in (|2.6| ) cancels. In particular, 



T 2i+1 (u,E) = -±—det(6 i jE6 1 (u + 2{£ + l-i)r J )+ 

OiW v (2.9) 



+ 6 i ^Mu+m+l-i)v) + k j+ Mu-2(i-m) 1 < t ! 2l + l 

is holomorphic at u — 0. 

The full family of operators commuting with <So is generated by Baxter's Q-operator 
Q(u). Moreover, the operators Q(u) commute with all the transfer matrices and among 
themselves: [T s (u), Q(u')] = 0, [Q(u), Q(u')] = 0. They obey the famous Baxter T-Q- 
relation 

9 x {u - 2£r])Q(u + 2?]) + 9 x (u + 2£rj)Q(u - 2?]) = Ti(u)Q(u) . (2.10) 
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We also recall the formula for T s (u, E) through eigenvalues Q(u) of the Q-operator [H], [l9| : 

Q(u+(s + l)r))Q(u-(s + l)r)) 



T s (u,E) 



U 2 p t 1 s 1 (u + (2£+l-s-2p)r 1 ) 



« Tlf-i e 1 (u+{2£+l + s- 2j - 2q)r}) 

X , V — \ — T 1 — — r , 1< S < 2£ . 2.11 

If s > 2£, there is no denominator in the prefactor. 

Let £ G Z + . In this case a commuting family, equivalent to the Q-operator, was 
explicitly constructed in [[L2|. Consider the operators 



A x = Y^ A k(x,X)e 

k=0 



(2k V -eri+\)d x 



(2.12) 



where A G C and 
A k (x, A) 



v ; [2£]! 



^ 1 0i(2x + 2(£-j)r l )e 1 (2X + 2(£-j» 

3=0 



#i(2.t + 2A + 2(£;-j» 



x g x (2z - 2(t-j)r,)6 1 (2X - 2(t-j) V ) 



3=0 



6x(2x + 2\ + 2(k+j-£)r)) 



(2.13) 



If = or k = £, the second (respectively, the first) product is absent. Here and below 
we use the "elliptic factorial" and "elliptic binomial" notation: 



[n]! = ni>1, \j) = ^jrj) , 

3=1 



(2.14) 



n 
rn 



n 



[m] ! [n — m] ! 



The main property of the operators ( |2.12| ) proved in [0 is their commutativity for all 
values of A: [A\, Ay] = 0. For generic A the chain of shifts in ( 2.12Q starts from — £i] + A 
and all the £+1 coefficients in ( J2.12 ) are non-zero. However, for X = It], I = £,£—1, . . . ,—£ 
only £ — \l\ + 1 of them are non-zero. (For example, Ae v = 1, A^-i) n = ([£]/[2£])So.) 

By difference operator in the next theorem we mean a finite sum J2k fk(x)e kridx with 
integer k. (So A\ are difference operators, in this sense, only if A = mq with integer m.) 



Theorem 2.1 ^71/ The ring of difference operators commuting with L = iS ^ s 9 en ' 

erated by L and A = A^ + xj v . 

It is convenient to introduce the following special notation: A = A^ + i) v , A = A-y+x)^, 
W = A— A. The role of the operator W (the very one entering eq. ( |1.4j )) for representations 
of the Sklyanin algebra is clarified in Sect. 5. 

We conclude this section by listing some properties of the operators Ax which will be 
useful in the sequel. 
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a) The Baxter T-Q-relation. It has been proved in [12] that the A x obey the following 
operator identity: 



S A X 



9i(2\-2£r)) 
0i(2X) 



A 



A+ry 



9 1 (2X + 2£r ] ) 
0i(2X) 



A 



\—r) 



(2.15) 



which allows us to identify A x with the Q-operator: A x = Q(2X) (cf. fl2.10|) , ( |2.4|) ). In 
other words, A x can be regarded as an operator solution to the Baxter relation. The second 
operator solution to the second order equation (|2.15|) is A_ x . Their wronskian Wr(A) is 
easily evaluated: 



Wr(A) = A x+rt A^ x - A\A- 



(A+r?) 



{[2£]\ 



,-1 



II 9 l (2X + 2j V )]W. (2.16) 



b) The symmetry x «-> A. Let F(x) be an arbitrary function. Since A^{x, A) = A^X, x), 
it is clear from ( |2.12| ) that the result of action of the Ax on the F(x) is symmetric under 
the interchange of x and A, i.e., 

A x (x, d x )F(x) = A X (X, d x )F(X) , (2.17) 

where A x (x, d x ) (respectively, A x (X,d x )) acts on the function of x (respectively, of A). 

c) Even and odd difference operators. Let S be the operator changing the sign of x: 
SF(x) = F(—x). It is clear from the definition that HA^S -1 = A_ x . We call difference 
operators O such that SOS -1 = O (respectively, SOS -1 = —O) even (respectively, odd) 
operators. It can be proved |T2| that for generic r\ any even difference operator commuting 
with L is a polynomial in L. In particular, A^ + A_k n for fceZ and A x A_ x for arbitrary 
A, rj G C are polynomial functions of L. 

d) Relations between the transfer matrices T s (u) and the difference operators A sr] , s G Z. 
From ( 2.13|) we immediately conclude that A_j V 



A ■ —I 



< 3 < 



so they are even 



operators. Then it follows from the above that the operators Aj V with integer — £ < j < £ 
are polynomial functions of S . So, similarly to (|2.8|) , we define polynomials Aj V (E) to be 
eigenvalues of the Aj V on their common eigenfunction \l/ such that = E^> . Comparing 
the fusion relation fl2.5|) for s <2£ with (|2.15| ), we identify 



A 



(E-s)ri 



[2£-s}\ 



T s (2£ V -(s-1)t] 



8 = 0,1, 



2£. 



(2.18) 



whence (|2.6| ) yields the determinant representation of the polynomials a v (E): 



A 



(i-s)r) 



' £ ' 




' 2£ ' 


s 




s 



det ES,. 



[2£+2- 



(here < s < 
polynomials: 



,J [£+1 

The Baxter equation ( p.15 ) gives the recurrence relation for these 



A 



(is-l)r)(E) = j^—^EA^E) + ^—rA 



(E) 



(2.20) 



with the initial conditions A er ,(E) = 1, A {e _ 1)r} (E) = ([£}/[2£})E. It is clear from ( ^20|) 
that 



.4 



(i-s)r,\ 



-E) 



■1YA 



(£-s)r. 



(E) 



< s < 



(2.21) 
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At last, we point out the relation 

T 2 e+i (0) = [2£] ! (A {e+1)v + A- {e+1)tl ) (2.22) 

which follows e.g. from ( |2.11|) . The difference operator in the right hand side is even. Its 
eigenvalue is given by the polynomial T 2 e+i(0, E) = lim n ^ T 2 £ + i(u, E) (see Q2.9|)). 



3 Double-Bloch eigenfunctions and explicit form of 
the spectral curve 

Let £ be a positive integer. For our current purpose it is more convenient to pass to the 
function 

i>(x) = *(x) ^J[9 1 (2x-2jr ] ) ) j . (3.1) 

Then the eigenvalue equation for the L acquires the form 

. 6 l (2x + 2£ri)e l (2x-2(t + l)rj) , . . 

which we also call the difference analogue of the Lame equation. In this form, the coefficient 
function is double-periodic with periods \ and |, so it is natural to look for solutions in the 
class of double- Bloch functions [[J, i.e., such that ip(x + |) = Biip(x), ip(x + ~t) = B T ip(x) 
with some constants B\, B T . 

Introduce the function 

Its monodromy properties in x are $(x + 1, Q = $(x, (), $(a; + r, Q = e 27r ^$(x, (), i.e., 
it is a double-Bloch function. The double-Bloch ansatz for the ip is 

1>(x) = £ s i(C, K, E)<b{2x - 2 3V , C) , (3.4) 

3=1 

where (, K parametrize the Bloch multipliers of the function vp(x): B\ = K*n, B T = 
K^e~ 2mC -. The coefficients Sj depend on the indicated parameters only. 

Substituting ( |3.4| ) into ( |3.2| ) and computing the residues at the points x = jrj, j = 
0, ■■■,£, we get £ + 1 linear equations 

i 

"£M ijSj = 0, i = 0,l,...,£, (3.5) 
j=i 

for £ unknowns Sj. Matrix elements My of this system are: 

J J 61,(20' + 1 )»))»! (2ji)) ^ 

^-i »i(C-2(j-»+lW 9,(2(i + W(2(i-<-l).)) , r 
+ * WO « I (2,)« 1 (20-i+l),) <** -*•'>■ (3 ' 6) 
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Here i = 0,1,..., £, j = 1,2,...,£. The overdetermined system ( |3.5|) has nontrivial 
solutions if and only if rank of the rectangular matrix My is less than £. By and 
M^ 1 ' we denote £ x £ matrices obtained from M by deleting the rows with % = and 
i = 1, respectively. Then the set of parameters (,K,E for which eq. (|3.5| ) has solutions 
of the form ( |3.4j ) is determined by the system of two equations: det = det = 0. 
So the three parameters are constrained by two equations. They can be written out in 
a particularly compact form in terms of the family of polynomials ( |2.19| ) and the elliptic 
"binomial coefficients" ( [2.14| ). Expanding the determinants with respect to the first row 
and taking into account ( |2.19j ), ( |2.15| ), we come to the following statement. 



Theorem 3.1 The difference Lame equation ( \3.2\ ) has double-Bloch solutions of the form 
( \3.4 ) if and only if the spectral parameters (, K, E obey the equations 



3=0 



e+i 



A jv (E) = 



Y / (-iy K ^e l (C - 2jrj)6 x {2(j - 1)77) 

j=0 



l+l 
J 



A u _ 1)v (E) = 0. 



(3.7) 



where Aj v (E) are polynomials of (£— \j\)-th degree explicitly given by ( \2.1Q ). They coincide 
with eigenvalues of the commuting operators Aj v introduced in t\2.12( ) on their common 
eigenfunction \I/ such that = E^> . 



The equations ( |3.7| ) define a Riemann surface T, which covers the complex plane. The 
monodromy properties of the ^-function (see Appendix A) make it clear that this surface 
is invariant under the transformation 



K 



Re inir, _ 



(3.8) 



The factor of the T over this transformation is an algebraic curve T, which is a ramified 
covering of the elliptic curve with periods 1, r. It is clear from ( |3.7| ), ( |2.21| ) that the curve 
admits the involution 

(C> K, E) i — > (C, —K, —E) , (3.9) 

so the spectrum is symmetric with respect to the reflection E — > —E. Another result 
of 0, which is not so easy to see from (|3.7|) , is that the curve T is at the same time a 
hyperelliptic curve. 



Theorem 3.2 j^j The curve T is a hyperelliptic curve of genus g = 2£. The hyperelliptic 
involution is given by 



(C,K,E) ^ {AN V -t,K-\E) 



N = -1(1+1). 



(3.10) 



The points P = ((,K,E) G T of the curve parametrize double-Bloch solutions ip(x) = 
ip(x, P) to eq. ($J\ ), and the solution ip(x, P) corresponding to each point P G T is unique 
up to a constant multiplier. 
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To compare with 
as follows: ( = z + 2Nr], K = k 



we note that the variables z, k used in that paper are related to (, K 
1 {z-2r j y 1 



. In terms of z, k the coefficients in ( |3.7| ) 

k-\ 



d 1 {z + 2r ] ) i 

become elliptic functions of z and the hyperelliptic involution is z — > —z, k 

The edges of bands ±Ei in ( |1.3| ) are values of the function E = E(P) at the fixed 
points of the hyperelliptic involution. As is clear from ( |3.10| ), the fixed points lie above 



the points ( = 2Nr) + u a , where u a are the half-periods: uj\ = 0, uo 2 
U4 = \t. The corresponding values of K are determined from 



2' 



UJ 3 



ha+r) 



Corollary 3.1 Let S a , a = I, ... ,4 be the set of common roots of the polynomial equations 



J2e a (2(N-j) v ) 
3=0 



A jv (E) = 



(3-11) 



J20 a (2(N-j)r ] )6 1 (2(j-m 
3=0 



J 



A (] _ 1)v (E) = 0. 



where N = \t{t + 1), and 6 a are Jacobi 6 -functions. Then the set of the edges of bands 
±Ei is the union of Ua=i £a and its image under the reflection E — > —E. 



To obtain a more detailed information from equations (|3.7| ), one can try to eliminate 
E and obtain a single equation connecting the two Bloch multipliers of the function ( |3.4| ) 
(parametrizing through ( and K). However, this is not easy to do directly. A possible way 
out relies on the following simple lemma. 



Lemma 3.1 Let be any solution to the equation 

0i (2a; - 2£ri) T . , 9U2X + 2£ri) T . , . . 

in the class of entire functions on the complex plane of the variable x, then 

V0v) = *(-3V), j = l,2,...,£. 



(3.12) 



(3.13) 



This assertion follows from the specific form of the coefficients of eq. ( |3.12| ). Indeed, putting 
x = in ( p. 12 ), we have ^(rj) = ^(—77). The proof can be completed by induction. At 
x = ±£r] one of the coefficients in the l.h.s. of ( |3.12| ) vanishes, so the chain of relations 
( |3.13j ) truncates at j — £. ■ 



Remark The conditions (|3.13|) resemble the " glueing conditions" for the Baker- Akhiezer 
function on rational curves with double points, where they are imposed on the \l/ with 
respect to its spectral parameter. However, contrary to that case, ( |3.13|) is imposed in the 
x-plane. 

Remarkably, the conditions ( |3.13| ) and the ansatz 



¥(a;) = K*'* [ I] 9 1 (2x - 2 3V ) ^ s m (K, C)$(2x - 2m V , () (3.14) 

U=l / m=l 
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for \1/ (equivalent to the ansatz ( |3.4|) for ip) with the same function $(x, z) given by ( j3.3|) 
allow one to find the relation between the Bloch multipliers even without explicit use of 
the difference Lame equation ( |3.12| ). Plugging ( |3.14j ) into ( |3.13| ), we obtain C equalities 
(for m= 1,2, ...,£): 

K™s m = {-lYK-™9 x {±m V ) ( n m ^2(m-j>) ) t^(m + n) V , c)s n . (3.15) 

This is a system of linear homogeneous equations for s n . It has nontrivial solutions if and 
only if its determinant is equal to zero, whence we obtain the equation for ( and K: 

det(K 2m 5 mn + G mn (0)^ =0, (3.16) 

V / l<m,n<Z 

where 

G mn (C) = (-l) £+1 [2m] ( fl [ f^)$(-2(m + n)ri,0. 

This equation defines a curve T e , which is the image of the spectral curve F under the 
projection r — > T e that takes ((, K, E) to (£, K). A more explicit description of the curve 
T e is given by the following proposition. 



Proposition 3.1 The equation of the spectral curve ( \3.l6j) can be represented in the f< 



form 



N 

E(-lYc}%)e l (C-4jr,)K« ir -fi = 0, (3.17) 

where N = \t{t + 1) and Cf\r]) are coefficients depending only on rj (and r) such that 
C?(v) = C%l j ( V ),Cl>%) = l. 

Proof. The following direct proof allows us to find the explicit form of the coefficients C- . 
To expand the determinant in powers of K, we make use of the identity 

f e l (x i +x j +o \ ^r 1 (c)^(c + 2Er = i^) A fffo-s,-) 

6 ^ B 1 (x i +x j ) nr=i^i(2x.) (H(.r, ■ 

which is a particular case of the formula for the elliptic Cauchy determinant. Let A be 
the set {1,2,..., £}. We use the following notation. For any subset J C A, A \ J is its 
complement, |J| is the number of its elements, and ||J|| = J2meJ m - Setting x n = —2i]n, 
we have: 



det 0~mn G mn 



l<m,n<e 

N 



J 

where 



E n n |?|±|M, (3.i, 

o y HU JCA,||j||=j fceJfc'eA\J °iK z K K K )V) 



k(J) = \J\£+±\J\(\J\-1). 



11 



Thus, the coefficient Cj reads 



JCA,||j||=i keJk'eA\j 
and the symmetry j iV — j is evident . 



Note that the sum in ( |3.19 ) runs over partitions of the number j into distinct parts not 



exceeding £. Examples are given in Appendix B. 

The meaning of ( |3.17| ) is the same as is explained after equations (|3.7|) : it defines a 
covering of the complex plane invariant under the map ( |3.8| ). This allows us to define the 
corresponding factor-curve which is precisely T e . Therefore, T e is a ramified covering of 
the elliptic curve with the modular parameter r. 

It is easily seen from ( |3 . 1 7|) that T e is invariant under the involution (£, K) i — > (4Nrj — 



C, K 1 ). Eq. (|3.12| ) allows one to express the function E through (, K by the following 
formula" 

_ 9 1 (U V ) g((l-l>7) _ ai-i t [1][2£] 

^(2^) *(^) [^-1] ' 1 ' J 

The coefficients s^, s^_i are given by the corresponding minors of the matrix K 2m S mn + 
G mn ((). It can be shown that E is invariant under the above involution of F e , so this 
involution coincides with ( |3.10| ). In terms of the ^/-function, the hyperelliptic involution 



takes ty(x) to ty(—x). Note also that eq. ( 3-17|) defines a singular curve. Indeed, the fixed 



points of the hyperelliptic involution are singular points of the curve T e , i.e., both the 



(- and ^-derivatives of the left hand side of (|3.17| ) at these points equal to zero. In the 
neighbourhoods of these points different sheets of the curve intersect. The function E 
takes different values on these sheets (which are obtained by resolving the indeterminacy 
in ( |3.2(J| )), so it resolves the singularities of the curve. 

Remark The function ty(x) is the common eigenfunction for all the commuting opera- 



tors A\ ( p.!2j ). Indeed, commutativity of L and Ax implies that ^\(x) = A\^(x) is an 
eigenfunction of L with the same eigenvalue E. By Theorem 3.2, \l/ is proportional to 
$?\(x) = g(\)^>(x). From the symmetry ( |2.17| ) and the normalization condition Ag v = 1 
we have 

Let us conclude this section by examining the behaviour of the spectral curve in the 
vicinity of its "infinite points", i.e., the points at which the function E has poles. From 
either ( |3.7| ) or (|3.17| ), ( |3.20| ) we conclude that there are two such points: oo + = (C — ► 



0, K — ■> oo, E —* oo) and oo_ = [C, — > ANrj, K — > 0, E — > oo). In the neighbourhood of 
oo± we have E = K ±x + o(-ft' ±1 ), while the leading terms of the function K are: 



K2= -W) n ir +0(lh c ^°- 
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The Baker- Akhiezer function 



^ba(x,P) 



p = (t,K,E)er. 



is easily seen to have the following asymptotics as P — > oo±: 

y BA (x,P)=K^ ±e (^(x)+0(K^) 



(3.21) 



with some functions ^(x), i.e., the pole divisor of the Baker- Akhiezer function is concen- 
trated at oo±. 



4 Explicit hyperelliptic realizations 

In this section we obtain the equation of the spectral curve of the difference Lame operator 
L, which has the explicit hyperelliptic form. This equation contains two variables: E and 
z. The latter is the eigenvalue of the operator = A, which commutes with L, 

on their common eigenfunction \P: = E^fr, Aty = z^. Recall that the operator 
A = A_(i4-i)^ commutes with both A and L. Let us write out the trivial identity A 2 — 
(A + A) A + AA = and act by both sides on the \1>. Taking into account (|2.22|) , we get 
z 2 - ([2£]\)- 1 T 2e+1 (0,E)+D 2e (E) = 0, or 



z + M9 = p]!)-i Ta+1 ( 0) E) 



(4.1) 



where D 2 e(E) is a polynomial of E of degree 2£ and T 2i+1 (0,E) = lim u ^ T 2 g + i(u, E) is 
the polynomial of E of degree 2£ + 1. They enjoy the properties D 2 i(—E) = D 2 e(E), 
T 2i+1 (0, -E) = -T 2i+1 (0,E). Recall that 



c(x)e^ +c+(x)e-* ) * 



where the coefficient functions are: 
0i(2x±eri) 



A = ^a 2k+l (x)e 

k=0 



(2k+l)rjd x 



(4.2) 



c±\x 



9i(2x) 



a 2k +i{x) = (-lY 



2£ + l 
£-k 



l-k-l 



2£ 



n 



B 1 (2x + 2(£-j)r ] ) 



k-l 

n 



B 1 {2x-2(£-j)r ] ) 



L 6 1 (2x + 2(£+k-j+l)rj) fJ- 9 1 (2x + 2(k+j+l)rj) 



(4.3) 

(see (£T2J), (^)). 

Introducing w = 2z - ([2£]!) _1 T 2 £ + i(0, E), we rewrite (|4.1|) in the customary hyperel- 
liptic form: 



w 



([2£]\)-\T 2e+1 (0,E)y - 4D 2e (E) 



2£ 



2£+l 



(E 2 



(4.4) 
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where P 2i+ i(E 2 ) = Yl^ 1 (E 2 — Ef). Note that w is the eigenvalue of the operator W = 
A -A. 

To find out the explicit form of D 2 e(E), we make use of the following simple argument: 
given a commuting pair [L, A] = of operators of finite order, the spectral problem for 
L is reduced to an eigenvalue problem for a finite matrix with elements depending on 
the eigenvalue of A (the "spectral parameter"). Specifically, let be their common 

eigenf unction. We set 

tyj = ^j(x) = #(x + jrj) . (4.5) 
The equation A^/ = z^/ allows us to express \l/o an d ^21+2 through \l/ 2 , . . . , ^2^+1: 



*o = z l ^a 2k+1 (x)^ 2k+ i 

k=0 



(4.6) 



-1 , nt J^a-2k-i(x + v) T 

*2i + 2 = za 2 ? +1 (x + rj)*! - E- hr~^^ 

k=1 a2e+i{x + r]) 



Now the spectral problem L\l/ = E^ can be rewritten as a homogeneous linear system 
for ^2^+1 ■ Equating its determinant to zero, we get a relation between z and E, 

which is the equation of our spectral curve. Its independence of the value of x follows, 
eventually, from commutativity of L and A. 

To be more precise, introduce the vector-function ty(x) with components ^j(x) (see 
Q4.5D ) , j — 1 , . . . , 2£ + 1. The eigenvalue equation for A can be rewritten in the form 

(4.7) 



•qj(x -rf) = A(x, z)^f(x) , 
where the matrix A(x, z) reads 

( z~ 1 a\{x) 2 _1 a 3 (x) ... 
1 ... 



A(x, z) 



z 1 a 2m (x) \ 















(4- 



V ... 1 / 

— * — * 

The scalar equation = E^ is then rewritten in the matrix form L(x, z)^{x) = E^(x), 
where 

L(x, z) = C + (x)A(x, z) + C_(ar)A _1 (x + T},z), (4.9) 

C± = diag [c ± (x + 77), c±(x + 2rj), . . . , c±(x + (2£ + l)r/)| . 

So, we have assigned to the scalar operators A and L the (2£ + 1) x (2£ + l)-matrices A 
and L respectively. 
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Lemma 4.1 The commutativity condition [L,A] = is equivalent to the Lax equation 



L(x - 77, z)A(x, z) = A(x, z)L(x, z) (4-10) 

for arbitrary values of z. 



Proof. This equality holds identically for the rows from the second to the last one. The 
first row gives the set of relations for coefficients of the L, A, which are equivalent to their 
commutativity. ■ 



The system L(x, z)^{x) = E^>(x) has nontrivial solutions if and only if E and z are 
connected by the equation of the spectral curve: 



det(Lio(x, z) — ESii) = . 

V yV ' ' lJ J Ki,j <2£+l 



(4-11) 



It is easy to see that this determinant does not depend on x. Indeed, denote the left hand 



side of (|4TTTD by f(x) 



It follows from ( |4.9| ) that L TOn (x+|, 



(x. z) . L mn 2 t 1 



whence fix) is a double-periodic function of the variable y = 2x with 
periods 1 and r and with finite number of possible poles. At the same time, ( |4. XCJ| ) implies 
that f(x) has one and the same value at infinite number of points x + mr], m e Z. Thus, 
f(x) = const. 

Extracting the ^-dependence of the determinant in ( |4.11| ), we can write 



det(L,ij(x } z) — Edij^j = z(—l) 



+ F + z- 1 G. 



where F and G are polynomials of E. According to the above argument, they do not 
depend on x. It is convenient to evaluate F at x — > — (£ + 1)77. This should be done with 
some care because some matrix elements are singular at this point. Using ( |2.9| ), we find: 



2£+l 



IJ 9 1 (2x + 2k V )\ T 2£+1 (2x + 2(e+l)r],-E) 



+ two "unwanted" determinants. 



At x -> -{I + 1)77 the first term yields (-l) m ([£]!)- 2 T 2m (0, E), and each of the two 
"unwanted" terms tends to zero. The simplest determinant representation for G is obtained 
at x = £rj when almost all elements of the first row are equal to zero. Skipping the details, 
we present the result. 

The equation of the curve has the form ( |4.1| ), where 

-1 

det(D ii ) 1 < ij < 2 ,. (4.12) 



D 2£ (E) = (-!)■ 



■[21 + 1] 



2£ 



15 



The (2i x 2£)-matrix reads 



( -E 



[21+3] 



[2] 

im 



-E #L 



[3] 








[6H] 



-E ^ 



[4] 



o o 






\ 





M 

[3£] 



V 



'/ : 



— £ 



[3m] 



[2£] 
[3£] 

—E 



(4.13) 



where the coefficients d^k entering the last line are: 



2k 



r 



k [* + 2k] 


' 3£ ' 




' 2£ + 1 " 




' 2£+k+l ' 


[k] 


£ 




£ + k 




k 



-1 



k — l. 



It is not difficult to see that the element Y>2i,2i-i is equal to 

[4^+1] [£+l}[2£] 2 



d 



21 



L 3£+l] ' 
The polynomials D 21(E) for small values of £ are written out in Appendix B. 

The size of the determinant in the equation defining the curve can be reduced for the 
price of more complicated matrix elements. Let us mention two examples. 

The first one is 2x2 determinant representation, which is in certain sense "dual" to 
( |4.10| ). The duality means that the roles of A- and L-matrices are interchanged: according 
to the same scheme, the spectral problem for A is reduced to an eigenvalue problem for a 
2x2- matrix with elements depending on the eigenvalue E of L. Specifically, set 



( o 



C(x,E) 



\ 



c + (x + rj) 
V c_(x + r/) 



and rewrite the equation L\& = E^ in the form 
Introduce the "monodromy matrix" 



Ec_ l (x + f]) 
^{x + rj) \ 



(4.14) 



C(x,E) 



%(x) 
ty(x + rj) 



A{x,E) = Y,[ a2k+ ' {x) 
k=o \ 





a 2 k+i{x + rj) 



Yl C(x + jri) 

2k>j>0 



where the arrow indicates the ordered product of matrices, and a2k+i (x) are the same as 

( ^(x) \ 

. Then the spectral problem Afy = is rewritten as A(x,E) I ^r^j^^ J = 



m 
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W(x + 77) 



and the Lax equation 



A{x + 77, £)£(x, J5) = £(x, £).4(a;, £) (4.15) 



holds true provided the operators L and A commute. Note that now A(x, E) plays the 
role of the Lax matrix (cf. ( f4.10|) ). The equation of the spectral curve reads 

det(A(x,E) - z) = 0. (4.16) 

Due to the Lax equation ( [4. 15[ ) it does not depend on x. 

The second one is I x I determinant representation for the factor-curve obtained by 
factorization over the reflection E — > —E. In the operator language, this amounts to 
finding a polynomial relation between the commuting operators Aw^n an d A^ + 2)rj- As it 
follows from (|2.15|) , their eigenvalues (e and £, respectively) are connected with E,z by 
the formulas 



- [24 [2* -l] v \t-iMJ' [i] V P + i] 

The equation of the curve can be derived in the same way as ([4.1 1[). 



5 Remarks on representations of the Sklyanin algebra 

In this section we make a few remarks relating the above material to representation theory 
of the Sklyanin algebra. 

Let denote the generators of the Sklyanin algebra realized by difference operators as 
in ( |2.2| ). To make connections between the commuting family Ax (|2.12j) and representations 
of the Sklyanin algebra explicit, we begin with a simple reformulation of the Novikov 
equation [L, V] = for coefficients of an operator V commuting with L = Sq . Suppose 
v(x,y) is any solution to the equation 

Vv(x,y) = 0, (5.1) 

where V = Sq (x, d x ) — Sq (y, d y ). (The operator Sq(x, d x ) acts to the variable x, etc.) 
Then the operator 

V = Y,V j {x)e?* B ' (5.2) 

with the coefficients 



Vj{x) =v(x,x+jrj) Y[6 1 (2x + 2(j -k)r])6 1 {2x + 2{j + k)rj) 

\k=l 




commutes with Sq . In general, the operators V given by this construction are of infinite 
order, i.e. the sum in ( |5.2| ) is infinite. For the family of commuting operators Ax Q2.12 ) 
the sum is finite. This corresponds to some very particular solutions to eq. ( |5.1|) . 

Recall that the Sklyanin algebra can be realized by certain difference operators 
in two variables acting on the invariant subspace of solutions to eq. (|5.1|) . Taking into 
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account the isomorphism between difference operators commuting with the operator 
and meromorphic functions on its spectral curve F with poles only at oo ± , we conclude 
that the Sklyanin algebra acts in the space of such functions on T. It would be very 
interesting to find the explicit form of this action. 

In the rest of this section we present some results on the role of the commuting operators 
A\ in representations of the Sklyanin algebra. First, we show that they have the same 
invariant subspace as the generators of the Sklyanin algebra ( |2.2| ). Second, the operator 
W = A — A is shown to "intertwine" representations of spins £ and — £ — 1. 



Proposition 5.1 Let £ be a positive integer. Then the operators A\ ( 2.12^ ) preserve the 
space 0j. 

Sketch of proof. Let F(x) G 0^. Then the monodromy properties of F(x) = (A\F)(x) 
are the same as thous of ^-functions of order 4£ (this is easily seen from (|2.13| )). Next, 
a further inspection of ( |2.12j ), ( |2.13| ) shows that the condition F(x) = F(—x) is enough 
for cancellation of all poles of F(x). Therefore, F(x) G 0^. It remains to prove that 
F(x) actually belongs to 0+ £ . Set f(x) = F(x)-F(-x), then f(x) = A x (X,d x )F(X) - 
A_ X (X, dx)F(X), where the x <-> A symmetry fl2.17| ) is used. Since f(x) G 0^, it is enough 
to verify the equality f(x) = in 4£ points x = mi] + u a , m = !,...,£, where u a are 
the half-periods. This is easy to do if to recall the operator identity A mr] = A- mr) for 
m = 1 i. m 



Corollary 5.1 For any X G C the operator A\ — A _\ annihilates the space Q\ v 
Indeed, for F{x) G 0^ the function [A\ — A_\)F(x) is simultaneously odd and even. 



Remark Recall the involution ( |3.9| ) that changes the sign of E. This involution takes 

I — . Z7TX 

the eigenfunction ty{x) (see p.l4j) to e i ty(x). Clearly, the operators A\ preserve the space 
e i 0J as well, so A\— A_\ annihilates the space 0^©e n Qj £ spanned by eigenfunctions 
of the difference Lame operator at the edges of bands (cf. H). 

To formulate the next proposition, it is convenient to slightly modify the operator 
W = A — A. Let us introduce the operator 



W 



-V. 



2£ 



<PT\x)W : 



(5.3) 



where (pe(x) = ]l,=o^i(2^ + 2(j — £)rf). The explicit formula for W can be written in the 
form that has sense not only for £ G Z + but also for I G Z + + ^: 



W 



2£+l 

E( 

k=0 



2£+l 
k 



9 1 (2x + 2(2£ - 2k + l)r]) 



nf=o +1 0i(2x + 2 3V ) nj =1 e,(2x - 2f V ) 



(5.4) 



Here and below the dependence of the W on £ is not indicated explicitly. The following 
proposition is proved by a straightforward verification using some identities for the 6- 
functions. 
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Proposition 5.2 For £ G |Z +; £/ie operator W "intertwines" representations of spin £ 
and of spin — (£ + 1): 

St t - i >W = WSi i >, o = 0,..., 3. (5.5) 
The same intertwining relation can be written for the quantum L-operator (|2.1| ): L^^^W = 

w\S e) . 

Remark In case of the algebra sl 2 the intertwining operator between representations of 
spins £ and —£—1 (realized by differential operators in x) is (d/dx) 2e+1 . It annihilates the 
linear space of polynomials of degree < 2£. 

Note that the operator W is not invertible. By Corollary 5.1, for £ G Z + W annihilates 
the space Qj e (Be *> 6^ (see the remark after Corollary 5.1). As is mentioned above, this is 
precisely the space spanned by eigenfunctions of the difference Lame operator at the edges 
of bands. So, in this way we obtain another proof of the result of @ : the eigenfunctions of 
L at the edges of bands span a (4£ + 2)-dimensional functional subspace, which is invariant 
for all Sklyanin's operators S a . The corresponding (4£ + 2)-dimensional representation of 
the Sklyanin algebra is the direct sum of two equivalent (2£ + 1) -dimensional irreducible 
representations. 
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Appendix A. Theta-functions 

We use the following definition of the Jacobi ^-functions: 

Oi{x\r) = £ exp [mr{k + ^) 2 + 2m(x + -){k + - 



fcez 



9 2 (x\r) = ex P (iriT(k + - f + 2nix(k + -)) 
fcez ^ 2 2 ' 

9^(x\r) = ^2 ex P (jfiTk 2 + 2nixkj , 

fcgZ 

O^xlr) = £ exp (-Kirk 2 + 2ni(x H — )k \ . 



fcez 



(Al) 



Throughout the paper we write 9 a (x\r) = 9 a (x). The frequently used transformation 
properties for shifts by (half) periods are: 

9 a (x ± 1) = (-l) 5 ^ +5 ^9 a (x) , 9 a (x ±r) = (-l) 5 ^ +s ^ e -^ 2mx 9 a (x) , (A2) 
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e 1 {x±l) = ±e 2 {x), 

e x {x ± f) = ±ie~ 1 4 mr ^ ix 9 4 (x) , (A3) 

0i(z ± ^) = ±e- l i niT ^ mx 9 3 (x) , 

By n we denote the space of ^-functions of order n, i.e., entire functions F(x), x G C, 
such that 

F(x + 1) = F(x) , F(z + r) = (-l) n e -™ T - 2 ™ a; F(x) . (A4) 

It is easy to see that dim9 n = n. Let F(x) G Q n , then has a multiplicative 

representation of the form 

n n 

F(x) = c Y[ 0\{x - Xi) , ^Xi = 0, 

i=i i=i 

where c is a constant. Imposing, in addition to (|A4]), the condition F(—x) = F(x), we 
define the space 0+ C 0„ of even ^-functions of order n, which plays the important role 
in representations of the Sklyanin algebra. If n is an even number, then dimG+ = + 1. 



Appendix B 



Here we explicitly write out equation ( 3.17|) for small values of £. The "eliptic number" 
notation [n] = 9i(2nrj) is used. 

£=1: 6 1 (()K 2 -6 1 ((-4r ] ) = 0, 

t = 2: 9 1 (C)K e - Be^C-^K 4 + ^((S^K 2 - O^-Ur]) = , 



[1] 



[1] 



[1][2] 



+ T7T? 1 ^(C-1M^ 4 - fflS^C^O^fT 2 + 0i(C - 24r?) = 



[3][5], 
[I] 2 



[1][2] 



4: x (C)K 



20 



0!(C-4r/)ir 18 4 



[1][2] 

/ [4][5][7] [4][5] 2 [6] 
V [1] 2 [2] [1][2] 2 [3] 

n [3][4][6][7] 



[I] 2 [2] 
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0i(C-12,)^ + fSii + [5]2[7] ^- — 



VW[2][3] [1] 



d^-lQ^K 12 - 



e 1 {c-2^ n )K 



10 



+ 



[1] 2 [2] 2 
/[5][6][7] , [5] 2 [7] 



+ 



V[l][2][3] [I] 3 
[3] [5] [6] 



e 1 {(-24r ] )K 8 - 



/ [4][5][7] [4][5] 2 [6] ' 
[1]*[2] + [1][2]'[3], 



0i(C-287/)AT 6 + 
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Let us present the explicit form of the polynomials T 2 e+i(0, E) and D 2 g(E) (see (|4.1|)) 
for £ = 1,2. 



[1][4] 


[2] 4 


[2] [3] 


[I] 3 [3] 




[2] 4 \ 


{ 3 [2] 





1 ; K U2] [i] 3 ; 

= [1] f[2][5] /[2] 2 [7] [4] 2 [5] [2][8]\ [4][5] /[7] [5] \ \ 

4{ } [3] I [3] [4] ^[3][4]^[l][2][6] + [4][6]; [1][6] VJ5] W/ J ' 
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